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ABSTRACT 


In this thesis, two structured receivers for a 
pulsed radar for detecting a Sworling 1 -target in the 
presence of a rotating random dipole modell^ clutter are 
considered. The clutter model incorporated- ic representative 
of chaff and vegetation clutter characters sed by spectrum 
with a slow secular component. The structured receivers 
considered are the conventional receiver and the mismatched 
receiver. The performances of these rocei VQrs 3-re evaluated 
for 4 and 16 hits/sean and they are compared with the 
performance obtainable with a white-noise optimum receiver 
(in the absence of clutter ). Their perf ormances, for the 
case of 4 hits/scan are also compared wi' till t^at of the 
optimum receiver which is difficult to i-Eiplement * study 

confirms that for targets with large doppl er > degradation 
in performance of the sub- optimum receiver is small. Bor 
a slowly changing clutter environment, ill® requirements of 
ah adaptive version of the mismatched recei 7cr are worked out 



CHAPTER-1 


INTRODUCTION : 

, / 

, The signal received by a Radar due to its own 
transmission may he considered to consist of two components 
one arising ‘from the scattering of electromagnetic waves from 
the desired target when present and the other due to returns 
from all other objects illuminated by the radar beam. The 
latter component is referred to as clutter and can he further 
classified as ground clutter, sea clutter, second time around 
clutter, precipitation clutter, angel clutter and so on. The 
characteristics of clutter depend not only on the relevant 
properties of the objects under consideration but also on the 
Radar system parameters such as frequency, polarisation and 
bandwidth etc. 

It is well known that the optimum receiver for the 
detection of a target in the presence of thermal noise alone, 
consists of a filter matched to the signal return from the 
target [1 ]. In, the presence of clutter, the ■ of 

the above white-noise optimum receiver-usually called the 
conventional receiver can be considerably degraded. But the 
design and implementation of receivers that are optimum for 
the detection of targets in the presence of clutter are 
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rendered difficult by two features that are typical of many 
engineering situations. The first of them concerns the absence 
of reliable statistical models for the clutter. The second 
relates to the practical difficulties that are encountered 
in finding the solution to the associated integral equation 
even when the statistical characteristics of the target and 
clutter are .available. Therefore in practice efforts are 
directed to the design of structured sub-optimum receivers 
with differing degrees of adaptation to and assumptions about 
clutter characteristics. 

For the design of both optimum and suboptimum receivers 
when the characteristics of the clutter spectrum are not known 
a, priori, they are quite often estimated on line. For 
estimating the clutter spectrum some appropriate parametric 
model is assumed and the model parameters are estimated on line. 

The clutter samples are modelled as an M parameter 
Autoregressive (AR) process and these M parameters are 
estimated by applying a discrete Kalman filtering algorithm 
to the message and observation model in Bowyer [ 26 ] . The 
maximum entropy methods are used to estimate the clutter 
spiectrum in Yiswesvariah [1411 • A detailed survey and list of 
references of this method are contained in [14-3* Wong*. Reed 
et. al. [ 23] developed a clutter model assuming "the clutter 
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to be due to a collection of rotating random dipoles and 
this model was further studied by Haykins [17] and the details 
of a practical receiver based on this clutter model is reported 
in [19]. However the ROC is not obtained for this receiver. 

In this thesis, we assume ¥ong’ s model for clutter 
and consider some receiver schemes for rejection of such 
clutter and compute and compare their performances in the 
Chapters to follow. 

ORGANIZATION OF THE THESIS'S 

After a brief review of the previous related work 
in Chapter 1, Chapter 2 discusses the clutter, target models 
used in this thesis. The optimum receiver problem is formulated 
and the various integral equation methods are briefly reviewed. 
In Chapter 3 the optimum receiver problem is considered ,. ■ 
assuming the transmitted signal to be rectangular pulse train. 

An expression for the performance index of the optimum 
receiver is obtained. In Chapter 4 some relatively easily 
implement able receives are considered. The receivers 
considered are the eonvuntional receiver, the mismatched 
receiver and the discrete resolution receiver. In each case, 
expressions for the performance index are derived. Further, 



for the mismatched receiver an adaptive scheme is studied 
in some detail. A quantization scheme along with a procedure 
for the clutter model parameters is considered. In Chapter 5 
the performances of the various receivers proposed are 
compared. The appropriate assumptions involved in the 
numerical computations are discussed. Chapter 6 provides a 
summary of the studies carried out and some suggestions for 
future work. A list of programmes used is given in the 
Appendix A. 
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CHAPTER- 2 

TARGET . CLUTTER MODEL AED PROBLEM FORMULATION 

Target detection performance of a radar receiver 
evidently depends on the electromagnetic scattering properties 
of the target and its environment. But he cause of the complexity 
and uncertainties involved in obtaining such information, 
radar receivers are often designed to maximize the probability 
of target detection in the presence of clutter and thermal 
noise for a given probability of false alarm that can be 
tolerated. In such an approach a statistical description of 
the radar echoes from the target and its environment is 
adequate for receiver design [ 5 ). In the following, we provide, 
to serve as a ready reference, a brief summary of the radar 
target and clutter models which are utilized in the design of 
receivers proposed in this thesis. Following this, wo give 
a preliminary account of the target detection problem and 
indicate the nature of the dependence of its solution on the 
statistical characteristics of the target and clutter. 

2.1 TARGET MODEL. 

A convenient and often used classification of radar 
echoes or equivalently of radar targets is provided by their 
fading and time (delay) and /or frequency (doppler) spread 
characteristics relative to the time duration and bandwidth 
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of the transmitted waveform. When the width of the received 
pulse returned by the target is approximately same as the 
transmitted pulse i.e-, when there is no spreading of the pulse' 
due to the target, the target can be assumed to be a point 
target.. Further if the fluctuation rate of the target is 
sufficiently small so that over a period of time equal to the 
pulsewidth of the transmitted signal* the target can be assumed 
to be stationary then the target can be modelled as a slowly 
fluctuating point target. The above conditions underwhich a 
target can be modelled as a slowly fluctuating point target 
can also be expressed in terms of the relative magnitude of 
the time duration, bandwidth of the transmitted pulse compared 
to the target characteristics [ 4 ]« 

For the transmitted signal normally used in a pulsed 
radar the aircraft can be modelled as a slowly fluctuating 
point target [6]. So ve shall consider some statistical models 
for the slowly fluctuating point target [7], [ 24 ]. 

Since the radar cross section (RCS) of a target is 
a measure of the proportion of the incident energy transmitted 
back to the radar, the targets can be modelled in terms of 
their RCS. The .statistical models for the radar cross 
section, for the slowly fluctuating point target L is given 

v 


in table 2.1 
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Table 2.1 

Statistical models for the RCS of slowly fluctuating 

Fluctuation Rate 


Target 

Probability density of RCS 

Model 

. . 

Swerling .1 

exp ( - % ) for x^O 

Swerling 2 

> * 

Swerling 3 

exp ) for x ^0 

< 

X 

. Swerling 4 I 

» r 

_ J 



Scan- Scan 

Pulse-Pulse 

Scan- Scan 
Pulse-Pulse 

Swerling 1 and 2 target corresponds to the case when the 
target consists of a large number of scatterers of equal 
echoing area. 

Swerling 3,4 corresponds the case when there is a large 
reflector in addition to a number of small reflectors. 

For Swerling 1,3 targets the RCS of the target can be assumed 
to be constant over a period of one scan but it varies from, 
one sc^h: to another. In Swerling 2,4 targets the RCS for 
successive pulses are independent. 

Ttfe shall assume a Swerling 1 target in' this thesis. 
Hence the RCS of the target remains constant over one scan 
period. If we denote the complex random multiplier which 
relates the transmitted signal to the received signal component 
due to the target as b then the envelope / b / is a Rayleigh 
random variable whose moments are 
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® / 3 - V" “^T" (2.1) 

B [/*b / 2 ] = 2 c£ * ( 2.2 ) 

2 

The value of usually includes the antenna gains and 

path losses in addition to the RCS of the target [3 J. 

2.3 CLUTTER MODEL > 

The clutter can be modelled as a doubly dispersive 
target for the transmitted signal whose time duration is 
usually small compared to the reciprocal of the bandwidth of 
the scattering function of the random scatterers and the 
bandwidth is large compared to the reciprocal of the length 
of the scattering area. Unlike the target return the interference 
or the clutter return is due to objects located at different 
range s . 

If we transmit a signal whose complex envelope is 
f (t) the reflected signal from the range interval (p,p+dp) 
is given by 

(t, ^ f Et (t-P- ) 7(t- u ) d P ( 2.3) 

where b (t,p. ) is the reflection process which is normally 
considered to be complex gaussian process.. The return from 
the entire range interval is 



9 


(t) = 


a 

/ y E t f(t-n)t(t- •J-,u)4(i 


(2.4) 


/-s«/ 


n^-(t) is considered normally to be a sample function from 
% zero mean complex gaussian process. It can be 
characterized, by the covariance f miction 

% (t,u.)=a [ J/3 t fJt-n)7(u-K 1 )S(t- £- >t i) _ 

C 11.1 .• 


'Efu- » ^l.)d p.d ^ ] 


(2.5) 


S tjR* (t *^ B £ ^ (*~ 2“ 

.. f(u- p^)dpdp^' ' (2.6) 


If we make two assumptions : 

1) The returns from different intervals are 
statistically independent, 

2) The return from each interval is a sample function 

of a stationary, zero mean, complex gaussian random process* 
Then we can write . , 


B lT( t , p)b( u, ) 3 =^q R ( t-u, p)6 ( ) 


The function (t, ) is a two variable function that 

depends on the reflective properties of the clutter. The 

0N*. 

Fourier transform of (v,p, )' is defined as the scattering 
function 


* 4 ^ 

(i.e.) )=/ Ut,|i ) (2-3) 
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Physically Sj^(f, p ) represents the spectrum of the 
process h (t, p ). Substituting (2*7) in (2.6) 

)= f(t— p ) Kppi(^—u> p )f ( n—p'Jdp 

~ aC .... (2.9) 

For different clutter K^Vjp ) 'will be different. 

SANGE INVARIANT SCATTERING- FUNCTION 


We want to consider the specific clutter for which 
the value pf the scattering function given by (2.8) at 
any frequency shift f is the same for all ranges. 

(i.e.) (f, n )is^ R (f) ....(2.10) 


Hence the scattering function has uniform doppler profile. 
We can treat the scattering function to be infinite in 
length if the uniform doppler profile extends beyond the 
range of a possible target. by 1 where 1 is the length of 
the target [ 5 3 . 

When the clutter is infinite in extent, the 


clutter return is a sample function of a stationary process. 
Using (2.10) in (2.9) 

% ( t ’ 11 )=St Ji ~ ( ^ ,x v : f : > 9 J f Sa* t_u 5 

a Et / f(t- P) f *(u- ^) % E (t-u) dp 
= Et K^t-u) mt. U ) f (u- ^ ) dM- ....(2.11) 


■■ 
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” B t K DR ( t ~ u ) / f ( t + T ) xf ( t ') dt' 

.... (2.12) 

-••'diC 


where t-p, = t'+v 

(2.13) 

a nd u- jx = t* 

(2*14) 


Prom (2.13) and ( 2 .14 ) "i£=t-u dt f =d \x 

But since the autocorrelation function of the transmitted 

signal R (t)« jf (t»+T) if (t‘ ) dt». ..... (2.15) 

, v • •• 

(2^12) can he written as 

*K n (*,u)» B t K^ R (t-u) *£(t~u) (2.16) 

We have derived the covariance function corres- 
ponding to the clutter with range invariant scattering 
function in detail since the clutter model that use in this 
thesis has this characteristics. 

2.4 THIS HAHDOH DIPOLE CLUTTER MODEL -* 

Wong et al. [23] developed a clutter model 
a g Burning the clutter signal to ho due to a collection of 
scattorers considered to he rotating random dipoles. The 
basis for the. above model can he described as follows. 

When an KF pulse is transmitted from a radar toward a 
cloud of random scatterers moving about and reflecting 
energy independently of ono another and if in addition 
the cloud is assumed to have an overall drift velocity, the 
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echoes returned to the radar will arrive at a rate which 
depends upon the local density of the cloud. If the effects 
of multiple scattering and scatterer rotation are neglected 
the echo signal from a particular scatterer can he regarded 
as a Doppler shifted replica of the transmitted waveform. 
However , any change in the orientation of the scatterers 
can cause variations in the returned echo power and phase. 
Hence the effect of scatterer rotational motion is included 
in calculating the echo return. 

The time varying correlation function of the echo 
signal is derived in terms of the characteristics of the 
transmitted waveform, polarization and the distribution of 
the scatterers. For the case of linear transmit-1 inear 
receive polarization and the transmitted signal consisting 
of a coherent pulse train, if the echo rate varies only 
slowly -with time and if the observation interval is long 
compared with signal duration ( so that the clutter process 
becomes stationary ) then the normalized autocorrelation 
function of the echo signal (clutter signal ) is given by [22l 

K ( t-u )= [ 1+i cos ave''' bT2 Je , ” d * i:2+;5CT R('c) (217) 

n c _ 

a = 4ix V or 

b = 8it 2 c£ 

o = 2* ^ 

d = 2w^cr 2 


where 
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In this model the scatterers are treated as dipoles instead 
of points with variable cross sections. The dipole rotation 
frequency and the doppler shift are assumed to be gaussian 
distributed. V Qr , a r are the mean, variance of dipole 
rotation frequency V^, are the mean, variance of doppler 
shift respectively. 

Though this model is derived assuming a 1 chaff 
cloud 1 this model is equally valid for returns from branches, 
loaves, grass etc. having a rotational motion under the 
effects of wind forces. 

Comparing (2.17) with (2.16) this particular clutter 

has a range invariant scattering function with 

2 . 2 

E Dr( t ) == cosavo”^ T ) (2.18) 

we shall be using this model for the clutter in this thesis. 

The advantage with this model is that it is enough to know 
just four parameters to know the clutter characteristics. 

2.5 TYPICAL VALUES FOR THE CLUTTER MODEL PARAMETERS. 

In order to obtain some typical values for the 
clutter model parameters we have to specify both the interpulse 
period Tp and the operating frequency f Q . The maximum values 
for, these model parameters are obtained by Haykins et< al 
[13], Barton [25]. Talhen the operating frequency is 1GHz the 
typical maximum values for the model parameters are j 
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dmax 


0.044 m / i 


T. 


V 


P 

0.1 


odmax T, 


m/ s 


0 

rmax 

Y 

ormax 


%£ 22 - hz 
P 


^222 hz 

-p 


For a PRF of 540 Hz the maximum values for a,b,c,d 

a = An V Qr Si 216 radians/sec 

b = 8it 2 a 2 Si 25720 ” 

c = 2-rc Y^ 524 1 » 

d = 11130 : (27t 2 cr d 2 ) »» 

* * 

2.6 PROBLEM FORMULATION . 

In the hypotheses - testing problem of detecting 
a radar target the hypotheses are : 

r (t) = B f Q (t) + n c (t) + < w (t) : H 1 

... . 

= n c (t) +'5f- (t) : H 0 

where r(t) is the complex envelope of the received waveform, 

f (t) is the dolayod, doppler shifted replica of tho transmitted 

signal, if is the complex gaussian random variable multiplier 

due to the target, ^(t) is the return duo to the clutter 

c 

and.'^(t) is the thermal noise assumed to be white and which 
includes the receiver noise. 
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The optimum receiver correlation waveform g (t) 
which maximizes the probability of target detection in the 
presence of clutter and thermal noise satisfies the integral 
equation [5]. 


T. 


r 

f s (t) 


- s 

Ti 




E^(t,u) 'g(u) du T ± *£t £ T f 




where n(t) = n (t) + w (t) 

t C 

and Eg<t,u) = B [ *n(t) n (u) ] 


K n (t,u) + !To'.6(t-u) 
c 


where (t,u) = B [ n' c (t) n c lu) ] 
c 

and E [ w (t) w’(u) ] = 11 o 6;(t-u) 


where Fo is the power spectral density of white noise. 

The decision is taken by computing the sufficient 

./•W* . 4*1+ . .A / 7 *2 

statistic f r (t) g (t)dt and then comparing 1 1Ll 

T. 

i 

with a threshold. ¥e want to consider the above hypotheses 
testing problem when the clutter covariance kernel is given 
by (2.17). This covariance Kernel is derived assuming the 
transmitted signal to be a pulse train. Ve ' shall consider the 
optimum receiver in chapter 3 and sub optimum receivers in 
chapter 4* Since we have to solve the integral equation in 
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oqn.(2.22) for arriving at the optimum receiver we shall 
briefly review the solution techniques for tho -integral 
equations. 

Substituting (2.25) in (2.22) we get 
T -f 

f g (t) = J K n (t,u) g(u) du + No g(t) /V* • (2.26) 

m C - 

x i 

This is a Fredholm equation of the second kind. 

2 • 7 REVIEW OF' METHODS FOR SOLUTION OF FREDHOLM INTEGRAL 
EQUATIONS^ . 

DEGEKEIfilTE ISIRTO,’' K3THCD . 

•* *• »***.. *. . « . * 

TjJhen the Kernel of the integral equation (2.22) 
is seperable into functions of t alone by u alone then the 
integral equation can be reduced to an algebraic equation 
[15]. Solving the algebraic equation* fee solution for the 
integral equation can bo obtained. 

STATE VARlAPir: 'APPROACH . 

Some integral equations can be solved by deriving 
a differential equation. Since tho differential equations 
are easy to solve the solution for the integral equation can 
be obtained by proper substitution [2]. For this to be 
possible the random process characterized by tho Kernel 
should permit state space representation. For the state space 
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representation to ‘be possible the Kernel has to be soperable 
in the variables t and u [!'). 

ORTHOGONAL SERIES EXPANSION METHOD. 

If the eigen values and eigen functions of the 
Kernel of the integral equation can bo found then tho solution 
to the integral equation can bo expressed in terms of these 
eigen values and oigon functions M. 

METHOD OP RESOLVENT KERNELS . 

The Resolvent Kernel of a Kernel is given by 

R K (t,T) = g. K in (t,T) (-l) n ~ whore K ±n (t,x) is tho 
n=l 

iterated Kernel given by K in (t,x) = J jr (t,s) K a _ m (s,T)ds 

T i 

and K 0 (t,v) = K n (t-x), In terms of the resolvent Kernel 

c 

the solution becomes [19] 

It 

g (t) = f B (t) t | ^ (t,x) f s (i) a-t 

1 

EOURIER TRANSFORM METHOD. 

In tho Integral equation if wo assume that the 
integration limits can bo extended to infinity tho Fourier 
Transform of both sides of equation (2,5) can bo taken. 


If we define 
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If the inverse Fourier transform of G(f) can be found then 
g(t) is known. 
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CHAPTER-5 

THE OPTIMUM RECEIVER 


In the previous chapter we had indicated that an 
integral equation has to he solved to obtain the optimum 
receiver correlation waveform. ¥e reviewed various methods 
for solving the integral equation. Further we also described 
the target and clutter models that are of interest to this 
thesis. To solve the integral equation, we need to specify 
the transmitted waveform. In this study we assume that the 
transmitted waveform shown in Fig. 5. 1(a) is a train of M 
number of equispaced rectangular pulse with pulsewidth T and 
interpulse period T^. 

3.1 AUTOCORRELATION OF THE TRANSMITTED WAVEFORM. 


The transmitted signal f (t) and the autocorrelation 

s 

function R(t) can be written as 
M-l 

f (t)=* 2T f(t-m^.) : _ .... (3.1) 

s m=0 

where f(t)= for O^t^T • ••• (3«2) 


= 0 otherwise 


R(t)= 


M-l 


k=0 


( 1 - 



) 


M-k 

M 


• •• « 


(3.5) 
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where p(T-kT p ) 

= 0 for 4T-kT p .f<^ T 
= j x-kQ}p| otherwise 


For M=2 

R(r) = X 

k-0 


2-k 

2 


/ _ P(\) 


= i- -4(1- 


P(t-ip) 


= 1- 


T 

LlL 

T 


T ” ) 

when j v j<C T 


• - ..... | t;_T | 

= £ (1- — — ) when (T~Tp(<T 


•••• (3.4) 


(3.3) 


5 . 2 THE KERNEL OF Tig, INTEGRAL EQUATION . 

The covariance function of the clutter signal can 
he obtained by substituting the value of R(t) from (3.3) and 
the value of K pp ( x ) from (2.18) into (2.16). We note that 
the covariance function of the clutter return is also the 
Kernel of the integral equation in 2.25. If f p is the doppler 
shift due to the target then (2,25) becomes 


jr fds-mljjei' V = 

m=0 



kL (t-u)g (u)du+Fog(t ) 
c 


...•(3*b) 


For simplicity we shall consider the single pulse case . 
(M=l) first. Then the Kernel of the integral equation (3.6) 
becomes 



J 
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2 

K n (l+JrCosa(t‘-u)e“^ > ^"' u ^)exp( 3 c(t-u) -d(t-u)^)x 

c ^ 

(1- %= ii, )^(t-u) .... (3.6(a)) 

where p(t-u)=l if (t-uj^T 

= 0 otherwise .... (3.6(.h)) 

Since (t-u) has quadratic terms in the exponent ((t-u; term) 

C # * | 

and has the factor (1 j~*-) it is not seperable in the 

variables t and u. Hence this is not a degenerate kernel. 

Since the Kernel is not seperable the state variable approach, 
and the orthogonal series expansion methods also fail. For 
tho Komel given by (3«6a) the iterated Kernels are difficult 
to find and hence tho resolvent Kernel method cannot be 
obtained. Wo must find the Fourier transform of the Kernel 
in order to apply the Fourier transform method to obtain the 
solution of the integral equation. 

3 . 3 CLUTTER POWER SPECTRAL DENSITY . 

To simplify tho notation we shall drop the tilde 

sign for indicating complex waveform. The clutter power 

spectral density S (f) is given by 

n c- 

-f-«C 

s (f) = Tk (•t)o“ 3w ' t aT .... (3.7) 

n c- *■' n c- 

-.00 

Putting t^u=T in (3.6a) and substituting in (3.7) 


we get o$ . 


S n (f)^J|"(l-+i co sai:e“ bT )e~<* T + 3 


CT 


P(t){ 1- -^)e- JW V 

=/ [e- a ' cZ+ j° T + e-( t>+a ) 1;2+ ^ a+0 ) T + i*eip(!' ' 


c 

oc 


06 


It I 




(b+d)Aj(c-ah)](l- -~- i ~) c | 3 (' 0 d'c .... ( 3 . 9 ) 


neglecting the normalization constant y and expanding. 


Since the term.' inside the square "bracket of (3.9) are similar 
excepting for the multiplying constants in the exponents let 
us concentrate on one term. Once we obtain the integral for 
the first term! the other integrals can be obtained by simple 
substitution. 

vXT 

Let S^ (f)= J e -dT 2 +jcT-3wT (1- i|i)p (T ) dT .... (3.10) 

C - — 9C 

For evaluating S^ (f) we shall consider two methods. 


3.4 FREQUENCY DOMAIN APPROACH TO EVALUATE! THE CLUTTER SPECTRUM . 


We can show' that 

O C 

JV tT 

— Off' 

Now let S’ 

X)R 


(l-'^)P(i)8‘ 3W1: dx=T( 


s in-rtf T\ 2 


oc 


rtfT 


r -s(f) 


.... (3.11) 


( f )= i 


-d-r +jcT. ^-jw-r 


dt 


— Oo 


« • • • 


(3.12) 




Then 




) =J e - dl:2 -j( w - 0 )' c aT 


<>c 


Completing the square in the exponent we get 
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%t< f ) = ez P [- 3 / exp[-d(T+j -gj2 ) 2 JdT .... (3.13) 

■ o»C 


«5C 


Since 


JV 


Q£3C j «• ir Ik 

dx= V rr~ 


-OC 


(3.13) "becomes 

p 

(v-c ) 

S^ R (f)= V-f- a” 4d 
*t*. os 

Now K' (x)=J (f )e*^ WT df 

n#% o 


.... (3.14) 


.... (3.15) 


Comparing (3*7) with (3.8) 


k; (t)^ +f”(i- if-')p(t) 

C 


.... (3.16) 


Since Z' (x) is the product of 2 time functions 
, n c 2 

p ( t)( 1 ^^) .and e“ di: + ^ CT (f) is the convolution of 


tho Fourier transform of c 

t i 

transform of (1- -jjj—) p(x) 


-d-r+5 


cx 


with the Fourier 


From (3.11) and (3.12) 

oC- 

S n J S'(x)££ R (f-x)dx 


— ac 

• <*5 


1 2,r SL - t*:gnr.q ) 


= j T( 0" V 4d dx ....(3.17) 


■ od " 
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Since the integrals j dx and f Q~ ax dx are evaluated 
using contour integrals we investigated whether contour 
integral can he used to evaluate (3. 17). ®qn. (3.17). can he 
rewritten as 


S n W= P* ii=S2|22Sl) f 0 -(f-z+fo) 2 k d .... (3 . 18) 

n c ^ (irxT) 2 d 


To evaluate (3.18) using contour integration we considered 


the complex function . 0 

TJT / _ . A • \ £ ' 


0(2 )=0 -K(*-f+f0)‘ 


2 2 
n z^T 


• • • • ( 3 » 19 ) 


whoro z is a complex variable* 0(,z ) has a single pole at 
z=0 and there is no other pole. Wo investigated evaluating 
the integral using circular contour first and then using 
rectangular contour. In both the cases the integral could 
not he evaluated since the integral is undefined over some 
segment of the contour { 9 [11]. 


3 . 5 INTEGRATION AND 


■inaaMafnaai 


I AT 10 


MHaiagaasaB 


FUNCTION. 


By integrating and differentiating the error 


function 


erf(t ) = yf J 


S e- 3 


.... (3.20) 
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oC 


and J cosbxe dx =V' 2 f ^ 2a (3.21) certain integrals 
e -t 2 

involving e can be evaluated [16]* For evaluating (3*18) 


ax 


we should consider the integrals 


r -k(x-f+fo) 2 r 

.. — dx an d I , 

*3T -s 


e 


*k(x-f+f o)‘ 


cosbx dx • . * . ( 3 . 27 ) 

X - x“ 

By differentiating and integrating (3*20) and (3.21) 
we cannot obtain the integrals of the form given in ( 3 * 27 ). 
Hence this method is not suitable, 

3.6 TIME DOMAIN APPROACH . 

We shall consider the method of obtaining S' (f) 

c 

from K* (x) directly instead of first finding S(f), 

V 

S^ R (f) and then convolving them together. 

Let us rewrite (3*10) as 

Bi (f>= f (1- ~t )P(^)e- D ' tZ -^ t+3oT dx(5- 28 > 

* o * ^ Oo* 

to avoid confusion of the differential dx with the exponent 
dx 2 in (3.10). 


Since 0(x)=O for (x\>T by (3.6b) 
0 2 
s; (f)= f (1«/I)e- DT 

II ^ 


.... ( 3*29 ) 


-T 


28 


T 


+ S (l-T/T)e“ I)T -3( w ~c)-r 
0 


T 


= 2R.P [ | (l-i:/T)e“ ;D ' T2 ~ ;5w,T d-rj .... (3.30) 

6 ' 


where w* = w-c 


Now (3-. 30) can "be rewritten as 


T 


S‘ (f)=2R.P[ i (1+ ^)e“' (T+jw ’/ 23))25) “' 4D 
0 

... w i2 


,2 


n 


dT 


T 


- -jjf J (T+jw/2D)e“* (T+;iw,/2D) D “ 4D dt] ..,.( 3 . 31 ) 
6 


T 


let I. 4 [(l +3W '/2DT)e- (T+3w ' /2I,)2l> - W ’ 2/41> ax .... (3.32) 
X 0 


and ( T+jw' /2D ) 2 D=z 2 


Then dirsidz/YD t=T ?> zsTyiJ+jw'/Pyi) 

. t= 0 ^ zsjw'^V^ 
TVD+jw'/2VD 2 

.'. Il=^g j 4 q“ z (l+jw'/^DTjdz 

dw*/2VD 


.... (3.33) 


... (3.34) 


ie(l+jw , /2M)[erf(TVl)+dw , /2Y^)-ei , f(dw , /2p) ] 

V - .... ( 3 . 35 ) 
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z -- 2 

where erf ( z )=2/)[% 1 e”" Z dz 

O'* 


(3.36) 


is the error function with complex arguement z, The values 
of erf(z)are tabulated in [ 13 ) . 


T 

let I 2 = j ( t+^W / 2D )e“^ T+ d w ' /2D)^D-w' ^/4! dT 

0 .... ( 3 . 37 ) 


and z=(t+;)W»/2D) 2 D 

Then dz=2D(T+jw , /2D)d-t 

2 

when t= 0 z=-w' /4D 

r=T z=D( T+jw* /2D)‘ 


.... (3.38) 


-w» 2 /4D 

V — X 

2 T 


D(T+ow»/2D)‘ 


-V 2 

■Ml 111 * 

4D 


1 f -z 
2D 


dz .... (3.39) 


w 
e 


,2 


D(T+jW/2D) 2 x 


asr [ e- 15 ) i2 ' ' = feTt ®' 151 _lw ' T -n 


-w 1 

4D 


.... (3.40) 


• S' (f )=2R.P[l 1 e“ w,2 / 4I>+I 2^ 
' n c- 1 
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= ^ E . P [ e rf ( IlfS + 3 wV 2 VS ) 34 - ( JjJA I ffi [ erf ( TVl >+ ff £)) 

+ fjj *" ^ l ) I m t arf ^ w '/ 21 fD )]- p E . P [ srf ( 3 w / 2 VB )] 

1 p — M : 

+ 2 DT *- Q co swT - 1 ] .... ( 3 . 41 ) 

Next lot us consider the case when there are two pulses in 
tho transmitted waveform . The autocorrelation function of 
the transmitted pulse train R ( t ) is given by 

R ( t )=1 - if hl<T 

' • f t-T f 

= i (1- ^ It-Tp! 4T .... (3.42) 

= 0 otherwise 

let S "( f )=/ R ( x) e 3o ^ 2 -^"dT 
n *. -00 

= ....( 4 . 43 ) 

—O0 

where R ( t ) is given hy ( 3 * 42 ). 

Substituting ( 3 . 42 ) in ( 3 . 43 ) wo got 

S * ' ( £ )= 2 R . P [ ! ( 1 - in ) a - 3 w -- i > x 2 ax H 

n c - 0 • 

0 - 3 w ' T - DT 2 aT 


I P+T 

J O- 


Hj.-Ti 
T ) 


Tp ' 


-T 


* * • • 


( 3 . 44 ) 
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1. 


p 


Let I = J (1— ?- + dt .... (3.45) 


T 

i P~T 
‘Pr. T 


Then I 3 = jjjE - ^ )e -D(T+j W '/2I)) 2 - W ’ 2 /4B ^ 


l p-t 

1 


(T+jw 


/2D)e- D ‘ 1:+ j w '/ 2D ) 2 - w,2 Al' a* 


(3.46) 


I 


P-T 


Ip ! 2 // T) 

= (1- Tjr -oW/2DI)e“ w /4D [erf(T p fD+jw'/2p)~erf 
(IpVD-TVP +jw'/2p)3- ^ p [e” T P D ~ dwtT P_ exp 
[~TpD-T^D+2TTpD+;jw 1 1-jw’ T p ] .... (3.47) 


Since the pulse width I <jtinterpulse period T p 
• T 

I,=(l- ip [erf(T T ,VB+iw'/2VB)-erf(T p Vli-TVB+3’' , /2Vl>)3 


.... (3.48) 


- m [i-« + 3w ’ t 3 


Similarly 

m p 

1 4= jf (1+ f“ )e~' D ' t ”^ W T df 

T P 
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(1+ +3W , /2DT)[erf(T p yD+TYD + 3wV2V^)-erf(r p V'D+3W»/2V^)3 

* » • » ( 3 * 49 ) 




2DT 


g-jw’T 


-lj 


*’• s ^’ (f )=2.A.R.P[I 5 3+2.A.R.P[I 4 3+S^ {f)+ ^ T [2coswT-2l coswT 3 


Tt\ sy 


(3.50) 


= 2A.R.P[erf (x+jy )}( 1- ^-)+ JyJ I m [erf(x+jy )] 

T 

“4(1+ ^)2.R.p[erf(x~Tyi)+jy)3- Ijerf (x+Tp+jy )] 

T. 


“4(1+ ^-)2.R,p[erf(x+jy))+ I m (erf(x+jy))- 

[cosw’Tp-cosw* (3} p -!T)) 

# * - 

+ f^[c°sw» (T p +T)~ c °ew , T p ]+S^ (f ) .... (3.51) 

c 

where x=T p VD> y=w'/2\fD, A=0(n)e“ w, ^ 4I) and 0(n)= 


Por S*'(f) n=l 

S ; a (f ), the spectrum for any arbitrary m can be found by 
c 

substituting T p by ml p wherever T p occurs in the expression 

for S*'(f) and replacing (f) by S^ m “ 1 (f). We observe that 
n c c c 
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the expressions for £ n (f) both for a single pulse and 

c 

multiple pulses are unwieldy. Hence we attempt to obtain 
an* approximate clutter spectrum. in the next section. 

3.7 APPROXIMATE CLUTTER SPECTRUM . 

Approximate expression for S (f) when the 

n c 

transmitted signal consists of a train of eqtiispaced 
rectangular pulses can be derived as follows ; The observation 
that the clutter is a slowly varying process [ 2 3 1, so 
that over a period of time T= pulsewidth of the transmitted 
signal pulses which, are of the order of p.sec it can be 
assumed to be constant simplifies the problem, 

let the number of pulses be M. The transmitted 
signal and the '.autocorrelation function of it are given in 
equations (3.1) to (3.5). We observe that R( t)=R(-'p) . 

When the received signal is heterodyned with an oscillator 
signal of frequency fo+fd whore fo,fd are the transmitted 
frequency, mean doppler frequency of the clutter return 
respo ctivoly 


,, , -bx 2 1 -(b+a)x 2 -jax 1 -(b+d)x 2 +jax' 

(t)=o + ~T G + 4 


then E 


DR 


(3.52) 


* • t • 
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a n& E dr ( x J-KppC -x ) 


*c 


S n 2,;R,P X R ( T ) E bR(' T ) e *”^ WTdT 


— *>6 
T 


M-l iT, 


- 2R 


t ziZz 1 

•p[ J Sjattld- «+ 57 [ J 

0 1=1 IT- 


z o>r( t ) 


P-T 


IT. 


(1+ $ ^ )o“ 3wt dT 


IT +T 
P +i 


1T t 


+ J W T)(1 " ? + -“Sr)®" 3 "* 

/ 

.... (3.53) 


IT, 


h T e shall assume that over a period of T sec. Kq R (t) can 
he assumed to he constant. 


M-l 1 T t 


Then S n (f)=2.R.P[E DR (0)/(l- dt+ lL 

C o 1=1 


J K dr ( lT p ) 


IT. 


IT iTp+I IT 

(1+ l - -jE).-!" d»+ J H^d^xx- f + yi) 


P-T 


it 


* ^ dTj (^) 3 


(3.54) 
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M-l 


2.E.P[ ^ (^JEjjjUT )I( - "| f - ) 2 cos2j.fl!r ....(3.55) 


i=0 


M-l 


= 2I( Si|2£L_)2 ^ [o-‘(iIp)‘ cos2 , flV 1 ,-OHdXlSp) 


1=0 


p + J 


[cos2*n:(f+f > )iT p +cos27i(f-f> )iT p ]] . **** ^ 5 ’ 56 ) 


where i t =B./2%i. 


3*8 VALIDITY OP APPROXIMATION . 

1X P+T 

¥e encountered integrals of the form J* R(v) 

iT 

K pp ( ^ )°""^ WT dv. ¥e stated earlier that Kj^t) cal. be 
assumed to ho constant over this interval. ¥e shall verify 
this statement, 

VJhen the PRP=300 T p =l/300 

Substituting in (2.27) 

c^max =13.2 c r max=V r max= 9 . 9 

a =12x9.9 ~ 120 
bCT-8000 

d= 2x10x169 - 3400 

K de ( T)= e - 3400T2 + \ ( 9 -lW00T 2 -ja V -llW + aa^ 
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T 


For simplicity wo shall consider only one term 
Lot 

For ir=iT p +T K£ R ( i^ p +T>e" 3400 ^ iT p +T ) 2 

= e -3400[(iT p ) 2 +2iTT p ] 

(VT^iTp) 


For T=iTp E^ E (iTp) = e ”3400(ilp) 


• K^i-Tp) 

' ‘ ZTJJ i*p+T ) 


| 3400(2iTTp) 


This ratio and honce the variation over a period of T seconds 
is maximum when i is maximum. For a transmitted signal 
consisting of 10 pulses imax =10. 


■ V 11 ?) I - ' 

“ KjjptiTp+T) I mar- - 


-4 

i 3400x2xl0xl0*' 6 x3xl0“ 5 = q 2x1 ° = 1 


— fi 1 

whore T= 10 sec. and T p = sec. 


Next let K dr (t)= e 


-11400v* 


^Dr( 4-^-p ) 
Then KppC iTp+'T) 


max 


0 6xlO“ 4 v x 


aT = 120x10 


^ , cosa( iT p +T ) 


r** 1 


cosaT p i 
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For K ( t \-, fi ~340£hr^ **-lI400ir^ .. . J _. 

-cor iLp^^T;-e -fe cosaT, the variation over 


'' ^se< 

< 

. d -= - T 


a period of T sec (10“ 6 sec) is negligible. 

Oj max V max 

or 

T 


o,max ■ v max ■ 0 max 

Next let us consider 0 ^= - 4 — ■■ , V = — ££ — ,0 = — £ 

or 2 r 2 


Then a=60 b=2xl0^ d=850 


K 


DR 


(<c)=e~ 85Ch: cos60re“’ 2850 ' T " 


let KJ E (t)= ,-W 

/ 

Ihen FVirTi) I - e 8 50(2ira p L e 2- 5^10-5 - 1 

“r\ t-i ( P 1- max 


DR 


let K- h (t) - e 


-2850-c 


:• max 
2 


* iT 3? ) 


=e 


2850(2iTT p ) =e 7.5xlO’ 


r5 




max 


Proceeding further in a similar fashion it can he verified 

that the approximation is valid for all possibles values of 

cr,,cr and ¥ . 

d r or 


3.9 THE OPTIMUM RECEIVER CORRELATION FUNCTION . 

We have derived the expression for S n (f). The 

0 

Fourier transform of the transmitted signal is given by 


r .( f ) 




M— 1 

■ J i Wn i. m 

i=0 


e 


-DwT p 


(3.57) 
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using (3.56) and (3.57) 


V f ) 

s ( f )= T^TTTfS o 


m/ sinitf T \ 


M— 1 

SI e -J wT p 

i=0 




... (3.58) 


If ■ we can find the inverse Fourier transform of (3.58) 
then the solution g(t) corresponding to the optimum receiver 
is known. But (3.58) is irrational and is the ratio of sine 
function and hence closed form inverse function g(t) cannot 
he found. 

3.10 PERFORMANCE EVALUATION OF THE OPTIMUM RECEIVER . 

Though we could not obtain the inverse Fourier 
transform of ( 3 . 58 ) to obtain the optimum correlation 
waveform g(t) we .can compute the performance of the optimum 
receiver that is possible if we had indeed used the optimum 
receiver correlation waveform g(t). We can compute the 
performance using 3*58. 

If the false-alarm probability is Pj, and the 
probability of detection is P^ then Pjr=(Pjj) 1+ ^ ....( 3 . 59 ) 

and the performance index & is given by [ 3 3 

‘ f 

= ® r J ? s (t) g* (t)dt .... (3.60) 
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whore is the average received energy. 
How (5.60) can he rewritten as 


T 

= B r J f f s (t) g* ( t )at = S r Jr s (f) 0»(f)df ....(3.61) 


T i 




using (5.58) in (3.61) 




t 

= E , F (f ) 

r J s ' Ho + S*(f) 


df = 1 


! 


IV*)| 

Ho+S (f) 
n c 


df 


If the doppler shift due to the target is f^ then 



IV f - f D>l 

Ho+S* ( f ) 
n c - 


2 

. df 


• • • • ( 3 . 62 ) 


Though closod form solution for the integral in (3.62) does 
not exist, A can he numerically evaluated. We shall 
consider this procedure in chapter 5. 
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CH&PTER-4 


THE COffTOTIOME RE GE ITER 


In the previous chapter we observed that the 

solution of the integral equation corresponding to the optimum 

receiver is difficult to find. However the performance of 

the optimum receiver can he evaluated. In light of these 

observations we want to consider in this chapter some 

structured sub optimum receivers which are relatively .easy 
t o implement . 

4.1 THE CONVENTIONAL RECEIVER. 

The conventional receiver is the white-nbise optimum 
matched filter. 1^ Tp,fp are the delay and doppler corresponding 
to the desired target the matched filter impulse response ‘.is 
given by 

M— 1 

f*U-t) = f*(-t*kT p )e +;Sw b t for ^=0 .... (4.1) 

• k=0 . 1 u 


M-l 

= f* ('t+kTp 4 Tp)e + ^ w D t for ^ 0 ....(4.2) 

k=0 


M-l 

when the transmitted signal is f^(t)= f(t-kTp). 

k=0 

The performance is degraded when the noise is coloured. The 
performance of the receiver is given by, 
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A 2 /Hl]- 3 [j|| 2 /H n ] 

E[(fi 2 /E 0 ] 

T f 

where g = J> r*(t)f s (t) dt 


» • • • (4.3) 


» • • • (4.4) 


In. the following, we derive an expression for£x . 


For the transmitted waveform given in 3.1 i.Q- 
oquispaced rectangular M pulse train 


T 


Tj 


2 


M-l 

= f r*(t) ^ f(t-mT p )o +;jw D t dt .... (4.5) 
'** m=0 




» fr*(t)f(t)^V dt+ J r*(t)f(t-T p )e +j V dt +,... 
J 0 T p 

(M~l)T p +T 

= J i»(t)e + 3V f(t-JFI TJdt .... (4.6) 

. (M-l)T p 


Eow as the doppler shift in frequency fp due to the target 
is small compared to the transmitted pulse width I the phase 
factor e*”^ w D^ can he assumed to he constant over a period 
of the pulsewidth T in eqn. 4.6. 
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Tp+T 


T . 

.1 = J r«(t)f(t)dt+c + jVr J r*(t)f(t-T p )dt+e + 3< ; M - 1 )VD 


T t 


(M-l)Ep+T 


X 


r (t)f(t-JPIa’ ) dt 


• • • » (4.7) 


(M-l )Tp 


substituting (3.2) in (4.7) 


T p +T 


£ = Jfr* (t ) yij? dt+e+^V? J* i*(t) ^ dt+e 3 ?I “ lT P W D s 

0 " T-r, 


(M-l)T p +T 
...jf r*(t). 
(M-l)T p 


fT 


dt 


.... (4.8) 


iTp+T 


Lot iT P^ J" r*( t )dt 


.... ( 4 . 9i ) 


iT, 


Then (£ = * * C M-1 

^ c l”* °M-1^ * 

i-e- Z &CV where C=[c 0 c p»«* c jy[-i3 


u=[ 


T 

T ^IxM 


1 

7t 

i 

W 


i 

w 


W& . . ( 4 . 9b ) 


Mil 
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««|2| 2 u^ T C ' T CU = U T C* T CU 

s [tff 3= U T E[C^G]U = U T kU .... (4.10) 


where k(i,j ) =E[c*c..3 

Now jTp+T iTp+T 

BEcJ c 3 3*E[ J X r* (tMuJdtduy^ 1 .... (4.11) 

j^p Up 

whore y = e^*” 1 W D T P 

Bqn. (4.11) can he rewritten as (using 2.21) 

jTp+T iTp+T 

E[4<!p=B[ / J [lf ( f i ’(t-iI E1 )Y* 1 +n*(t)+w*(t)3 

[Bf(u-jTp)y^+n 0 (u)+w('u)]atduY : ' _:1 '] . .... (4.12) 

= [Jj [B itB|V‘(t-lI p )Y ,iri f(«-J®p)]+B[B > i ,, ''(t-iI p )Y* 1 ii c (u)] 

+ E[^ /‘(■t-iT p )Y 1 ' :L w(u) 1 +E[Bf(u- 3 Tp)Y‘ l n^(t ))+E[n*( t )n 0 (u)) 


+ B[n*(t)w(u)1+E[w*(t)Y j f(u-3Tp)B)+E[w*(t)il 0 (u)) 


+ E[v' r (t)w(u)] dtduj V’" 1 


(4.13) 


• • * • 
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If we assume that the clutter return and receiver white noise 
are uncorrelatod and assume that 


E[w*(t)w(u) ]= No6(t-u) 
B[n*(t)n^(.u)3«K^ (t,u) 


• • • • ( 4*14 ) 


Then (4.13) he comes 

E [ c^Cj ] =f [Xa ( t , u )+H 06 ( t -u )+f* ( t-iT p ) f ( u-j T p )f ‘ 3+1 
iB [jB| 2 ]dtdu .... (4.15) 

No6(t~u)dtdu =No6 (i-j)T .... (4.16a) 


j Tp+T iT p +T 


How j J 

3 T p iT p 


j T p+ T iT p +T 

and J | fnt-iT p )f(u-.jT p )Y ;5+ - i EC/B| 2 3dtdu 

jTp ±T P 
jT^+T iT^+T 

= X f (u-3T p ) J f*(t-iT p )f 3+i dtdu E[^Bi 2 3 ,...(4.16h) 


3*x 


iTr 


= i- . 4r .T.T.y 3 ~ i E[|Bl 2 3= Y" 3+a ®B[iB* 2 3 C4.17) 

y i V T 



45 


Hq-w 


3T p +T 

I 




• • • » ( 4 . 18 ) 


where R(t-u) is given "by (3.3) and is reproduced here for 
convenience 

M-l 

E(t-u) = 21 jpU-' 1 ? 0 ) .... (4.19) 

k=0 

■whore P(« k )= h-kT r | for 
= 0 otherwise 


. Since K j)^(" fc “ u ) is the lowpass clutter process is 
slowly varying and hence over a short period of time 
T=pulsewidth it can be assumed to be constant. We verified 
this in the last' chapter. 

Then 

jT p +T iT p +T jTp+T iT p +T 

f f Ejj R (t-u)R(t-u)dtdu = r w iTp- T ) f R( t -u )dt du 

- W _ - 


3T p iTp 




iT T 


jTp+T iTp+T 

a« Da ( i^p -3 Tp ) f J R( t-u )dtdu 


OTp iT T 


• » • • 


( 4 . 20 ) 
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let A(k) = |~£ (4.21) 

f' /** 

To evaluate J_J P( t-u)dtdu (i.e.)(4.20) we shall consider 
the various values of t,u and the corresponding values of 
k for which J and hence we shall find H(t-u). 

The values are given in table 4.1 


k= -[(t-u)modT ]+(t-u) 


•••• (4.22) 


No 


u 


t-u 


te(t-u) 


x-kT, 


iT P 

iT p HT 


iTp+T 


4 

5 

6 

7 

8 

9 

10 


iT p 

iT A 

it p +|t 

iT p +T 


IT • 

P 

ilp+T/2 


iT p +T 


' JTp 


• ' I I 


t 1 


«vfi 

t I 


I I 

I » 


W-3)V , 

|A(i-j)[l^| 3 


,01-3 


(i-3)T p +T 


r i 


I t 


|A( i-j ) [ 1-Tjr 3 


0 


T 

2 


T 


T 


(i~j)T p -± 


(i-i) 




T 


( i-j )T p +^ 
(i-3)T p -t|l 


I f 


t t 


t « 


A( i-3)[l-fp] 

A(i-j)[l- $ 3 
A(i-j )(l~ |jp] 

Ki-j)[i-§ |3 


T 

5 


T 

5 


^A\ < 


T 


(i-jl) A( i— j )[1- | T 3 


! I 


t ! 


(i-j)T p 


T 




I f 


» f 


A(i-3 )[1- § 3 

A(i-j)[l- fj] 


~ T 

3 

"!£“ 

5 


.0 


T/2 


TABtt + *l 
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Prom ‘1,2,3 entries in table 4.1 we infer that 


iTp+T 


T 


J R(t~jTp)dt =A(i-j)X (1- m ) dx 

jr ■ o 


iT. 


.... (4.23) 


Prom the entries numbered 4»5,6,7 we infer that 

f* R(t— jTp** ^)dt = 1( J (1+ iy)dx+ ^ ( 1~ ^)dx •••• (4*24) 

* - 2 m ' ■ v\ 


LT +T 


LT. 


-r 


o 


Prom the entries numbered 8,9>10 we get 
LTp+T ® 

jR(t-jT p ~ f )dt = A( i-j )[ J (1+ $) dx+ <J (!- f) dx 3 

i-Tp ' "2 O'. 

.... (4.25) 


combining (4.23)> (4.24) and (4.25) we get 

iTp+T 0 y 1 

f R( t~jT p -y)dt=A( i-j )[...J (1+ 7f) dx +.,F (1- ^) dx 1 ..*.(4.26) 

wi i. . 0 

-(T-y) 


iT. 


,jT p +T ±T +T S iT,+T 

. j ,:i R(t-u)dtdu = Jf J R( t -3Tp-y) d ' fcd y 


3V 

T 


0 iT. 


o y 

= f A( i-j ) [ J (Hf)dx + J (l-f) dx 3 d y .... (4.27) 

"0 ' -(T-y) 0 


substituting (4.28) in (4.20) we get 

jlp+T iTpT 

J /k (t.ujdta^Kj^UTp-jipK ft-jft- l L ] | 

iJL-rv v- 


T 


3T. 


.... (4.29) 


Substituting (4.16a), (4.17) and (4.29) in (4.15) we get 
E[c£ <> 4 ] =Ho6(i-j)T+ | T[ (i-j)T r Y (l_:i) +TE[»Bli 2 ] 


• • • . ( 4*30) 


let *L be a matrix whose elements are 

k l(i» 3 ) =E f c i c j^ 


./• E[ 2 /H 1 ) - u kpTJ 
Similarly B[fCl 2 /Hj=U T 'k U 


.... ( 4 . 31 ) 
.... ( 4 . 32 ) 
.... ( 4 . 33 ) 


where k Q is an MxM matrix whose (i,j)th elements are given by 

• ' * 0 U»3> X£U.3H*E B 2 ] (4.34) 

* • U T K p XJ 

« * & ~ !3E 


.... ( 4 . 35 ) 


U A K 0 U 


where K 2 is an MxM matrix whose (i,3) l* 19 element is given by 


*- 


£ 2 ( 1 , 3 )= TE[ > Bl _3 » 1» 3 = 0,1,2. . .M-l 


• • * • ( 4 * 35& ) 
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The Vector U is defined in eqn. 4.9. The expression for 
the performance index 4, derived in [ 14 3 also has the same 
form as (4.35) excepting that the matrices are d<Tf ' 1 * 1 ’ i ad 

and obtained using different statistics. 

AH' EXAMPLE 


The computation of the performance index can be 
made clear by considering a simple example, let the total 
number of transmitted pulses be three. Then the elements 
of the relevant matrices are 


K 0 (0,0)=NoT+ | T^O) 

K q ( 0,1)= |T(2=i)K DR (l)Y 

K q ( 2,2)= HoT+ | TK DR (0) 
U T K 2 U=[^r |f p) 


let E[Ib! 2 ]=1 

_1 = I T W 1 ) 


— - 


1 

T T T 


w 

T T T 


1 

Vt 

T T T 


i 

yr 


f £ 


ss 9 36 ) 
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U T E 0 U-[pr ^ p 


3 


[ HoT+ fTK^O) |-TK dr (1)y" 1 


9 

2 
9 


9 X2V DR' 

I TE DR (2 ^ 


FoT+ 0 


2 9 TK Er( 1 >Y 


. where y=e^~'^ w I)'^P 


^ K DR^ 2 

9 tk Dr( 1 )Y 

NoT+I^CO) 


1 

75 


7? 



=3Fo+ | 1^(0)+ ^■R DR (l)cosw ]D T p + | E DR (2)cosw D (2T p ) ....(4.57) 
substituting (4.36) and (4.57) into (4.35) we get . 


Prom (4.37) we observe that & is minimum when w p =0. Hence 
maximum degradation in performance occurs when the target 
doppler shift =0. As the doppler shift increases 6 also 
increases. In the next chapter we compare the performance of 
the conventional receiver with that of the optimum receiver. 


4.2 MISMATCHED RECEIVER . 

The conventional receiver is simple to implement. 

But the performance of the conventional receiver deteriorates 
as the clutter/white noise ratio increases and. when the 
relative velocity between the target and clutter is zero. The 
mismatched receiver is the receiver which retains the simplicity 
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of the conventional receiver but performs better than the 
conventional receiver when choson approj^rihtoly. 

. For a transmitted signal 

M-l 

of the form f^(t)= 2T f(t-KT p ) the mismatched receiver uses 

k=o - . * 

a correlation waveform consisting of constant amplitude 
(rectangular) uniformly spaced pulses with different weights 

V«S. 

* 

If the doppler shift due to the target is . f ^ the 

receiver correlation waveform > 

M-l • 

V(t)=2~ v,f( t-kTpJe^V .... (4.58a) 

k=o * - r 


The optimum v k *s which maximizes the performance of the receiver 
will be found in the next sections. 


4.3 PERFORMANCE INDEX OF THE MISMATCHED RECEIVER . 
The mismatched receiver computes 


£ m =/r(t) ^(t)dt .... (4.38b) 

and compares 2 with a threshold. v(t) is the arbitrary 
function that we want to choose [53 

The performance of the receiver is specified by 


j^-m where 


using (2.21) and (2.23) in (4.38b) 


• • • • (4. 39 ) 


1 t 


■ Nf-UH 
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n 

E Kmi 2 / H i^ =:E ^ ^ J/f Bf s ( t 5+3i(t)]v*(t)dt ) 2 ] 

TV 

” JjV u )v* ( t M u )d tdu+ j j E [Bf B ( t )n* ( u ) J v*( t )v( u )dtdu 

iJjE[B*f* (u)n* (t ) J v* ( t )v( u )d tdu+lj^ E[ l B | 2 ] f g ( t )f u 1 
v*(t)v(u)dtdu 


If the signal and noise arc uncorrelatod 

Then 

E tieJ/ H i3 *0 V t > > U ) T *( * ) v ( u)dtdu+fjf s ( t )E[» B » 2 3 f*( u)v( t )v(u)dtdu 

* • • • ( 4*40) 

Similarly |£ m |?/H 0 *(t)v(u)K a (t,u)dtdu ....(4*41) 


Substituting (4*40) and (4.41) in (4*39) wo get 
B[ |Bi 2 *f f ( t )v*( t )dt { 2 ] 


A as 

m 


J| v*(t)v(u)yt,u)dtdu 


* • * . ( 4 • 42 ) 


For tho uniform dopplor profile case 

K n (t f u)» J K^jj(t-u,)f g ( t-(i)f s (u-^)d }4 +No 6 (t~u) 

* K DR( t - u )J f s { t~ p)f e ( u- p)d p +Ho6( t-u) .... (4.43) 


- «4. 


Substituting (4.43) in (4.42) 

_ B[|Bl 2 ]|Jf B (t)v»(t)dt{ 2 

( t M U )^ H ( t-u)f B ( t- M)f s ( U- tfdtaua n 

(4.44) 


# ♦ • * 



53 


.4 MAXIMIZATION OF THE PERFORMANCE INDEX FOR MISMATCHED RECEIVER 
For maximizing iJk^wo can maximize IT' g (t )V* (t )dt| ^ 

2 e ping the denominator of constant. 

.e Maximize tfi (t)v*(t)dt> ^ subject to the constraint 

//v*(t)v(u)E n (t,u)dtdu =C* .... (4.45) 

sing the Lagrange multiplier yVwe shall form the function L 
aich has to be maximized without any constraint. 

.0 L=//f s (-t)-s»(t)at I 2 + A[//v *( t )v( u)K^( t ,u)dtdu-C3 

» . • • ( 4.46 ) 

ow (4.49) can be rewritten as 

-//V t)f g (u)v* (t )v(u)dtdu+^ [C-^Jjfv (t )v(u)X n (t ,u)dtdu 3 

.... (4.47) 

e want to maximize L by choosing v(t). Let v Q (t) be the 
unction which maximizes L w.r.t. the integration w.r.t the 
ariable t. Let us consider the effect of variation of y G (t) 

y€h(t). With 

" * * • ’ 

v(t) t=v Q ( t )+£h(t ) .... (4.48) 

et the value of L be denoted as I^* 
hen 

% s ( t )f B ( u) [v o ( t )+e h» ( t ) ] v(u )dtdu 

+ [C^ t fj[v 0 (t)+^h # (t)3K n (t,u)v(u)dtdu *... (4.49) 
f v^(t ) is indeed the optimum v(t) which maximizes 
. w.rlt t then the variational term in should be zero [ 10 ] 
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'-//v ^ ) v ( u)dudt- (t )K^( t ,u)v(u)dtdu =0 

.... (4.50) 

;i.e)^h*(t)£[ f s (t)f^(u)v(u)- E n (t,u)v(u)]dtdu=0 

.... (4.51) 

Cf this is to be true for any arbitrary h*(t) then it 
if ill be possible only when 

f s^ t ^X f s^ U ^ U ^ du ’" / ^T V fc,u ^ u)du =0 

i.e f s (t) Cl =^/v t ,u)v(u)du 

*** f s (t) * t , u )v( u )du=^/ K n ( t ,u)v(u)du ....(4.52) 

fhoro CjS* Jf*(u)v(u)du ^i= ^ .... (4.55) 

Similarly it can be shown that the optimum v(u) which maximizes 
j w.r.t integration w.r.t u also satisfies the same equation. 

\ I is maximum when v(t) and v(u) satisfies the equation (4.53). 

When we maximize with no contraint on v(t) the 
>ptimum v( t ) is the same as the optimum correlation function 
j(t) that wo shall use for a coloured noise. However when 
re restrict the v(t) to be of the form given in equation 
4. 38a). the solution of 4*53 for this v(.t) yields the optimum 
r alues for the weights. 
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4*5 DETERMINATION OF VEIGHT VECTOR. 

The correlation waveform v(t) corresponding to a 
mismatched receiver given by (4.38a) can bo rewritten in vector 

form as . 


v (^)- ? 2*** V M-1^ 

whoro y=Q~^ W D T P 

Prom (4.53) tho optimum values of 
v i l=o»l.* M-l can bo found as 

follows : 


f(t) 

f(t-T p ) Y 


f(t-mTp)Y 


M-l 


substituting (4.38a) and (4.1) in (4.53) 


M-l 

f{t-KC )e 3 V« ( ,X V t,u ^ v k f (^ kT p) e3Wl>U du 

k=o * _y k=o 


T, 


M-l 


T. 

l 


Tp+T 


*= j K n (t,u)v 0 f(u)du+ J K n (t,u)v 1 f(u-Tp)e3 w D T P u du + ... 

0 Tp 

( M-l )Tp+T 

+ ' J K n (t,u)V M _ 1 l( U -fPlI p )e :iw D (M - 1)T P u du ....(4.54) 

(M-l>T p 

«= ^ [G 0 (t) C^t) 3 


v 0 

V, 

V 0 Y ; 


W 


M-l 


. (4*55) 
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where Y=©‘ }W D T P 
iTp+T 

Ci(t )= Jf K n (‘ t f' u ) du 

iTp 

Let t= ^ Tp 4 £ ~ 0,1. • M-l 


. . .. (4.56) 


iTp+T 


iTp+T 


'.^(.^Tp) = j K n (t,u)du = j K (t,u)du+No6(t-u) 
±T iTp c- 


lT p 1J, c - 

"P P 

iTpifT 

•fx* (j?Tp-u)du +No6 (i-i) 
c 

iTp+T 


(|,!Ep-iTp) J R(^Tp-u )du+No6( i-£) 


iT, 


(* * KjjpCt-u) is slowly varying ) 
iTp+T , T 

How J~ R(£Tp-u)du=A(i-| r ) J (l- |)dx= ^ A(i-£) 


iT. 


M-( i- £ ) 
= ' ' 2 K ~ 


.... (4.57) 


Substituting (4.57) in .(4.56) we get 

c i (^T p )=E dr ( £-i )T p . +Ho6( i- i ) . . . . (4.58) 

- f(0)= ^ [C^(0) Cp(0) .... 0)3[ 

T 

However £(0)= yjj? 


V 0 

V 


W 


M-l 


.(4.56) 


f* • • • 


(4.59) 



57 


Similarly 

f(^p)= [C^(Tp) Cp(Tp) . C^__p(Tp)3 


and 


\y 


V M-1 y 


M- 3 ] 


( 4 . 60 ) 


,M-1- 


f(H=lTp) =^~* [C 0 (M~lTp)C 1 (MIlTp) T p )l V 

• • • • ( 4 . 61 ) 

From ( 4 . 59 )* ( 4 . 60 ) and ( 4 . 61 ) we infer the following relation 


1 


0 0 (0). . . 

Y 

Y 2 

• 

= 

• 0 0 (lp). . 

♦ 

•- 

+■ 

• __ 

y M-l 


♦ 

0 0 (H-lT p ; 

— — 1 

* 

i 


• • • • 


,\ 


W°> 

wv 


■fcl 


0 


vlY 

Y 0 y‘ 


Vi Y 


M-l 


“I ^ 

Y 0 

2 


H O 

r 

Y 


7 n 

o 


2 

• 

• 

o 


_^M-1 


* • • .( 4 . 62 ) 


....( 4 . 65 ) 


i.e F = £_ T* V where p = 


1 "I 


1 

Y 0 


• 

and Y ^ 

-* 

• 

0 . 

« 

M-l 

L° y i 

j 

■ ♦ 

v 

JM-lJ 

-IjT 

m • * • ( 4*64 } 
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V= r "g* P . . V= f'C'p .... (4.64a) 

The weight vector for the mismatched receiver has 
been obtained hy Rummlei[2G] [2l] , Spaf ford et al [ 2,2 ] using 
different methods f^>r maximation. The final result obtained 
is the same as that we have obtained here. 


AN EXAMPLE 

Let us consider a transmitted waveform consisting of 3 
rectangular pulses 


Then 




— — 


' 

10 0 

\ 

i ' 

O 

O 


Y 

0 0 Y 2 


! Y 2 

_ 


._j 


... (4.65) 




3(0,0)= C 0 (0) = Ho+ —ft — -08(2,2) = 3(1,1) 
E(0,l)=0 0 (T p )= KpgtTpH | - 3(1,0) 

^o,z)=e 0 (.2%)*K BB i z ' l! 2) ' i • \ =13(2,0) 

3(2,1)- Sjii’.j.Itl,!) . 


__ W 0) k dr (t p ) K m^hl 

ITO+ 1 7T 1 *— 3? 


... (4.66) 


r* ~ 

v o 

V 1 

rnimm 

10 0 i 

0 Y^O 


V 


0 0 Y~ 2 

i 


) l u:r - V- 





Ed R ( t p ) 

5 


TT k dr( °) 

N o 1 — 


K BR^ T P l % R 

g —y 


No+ 


r 

3' 

K DR^ 0 ^ 


-1 
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1 0 


0 y 
0 0 


-1 


— 






~ “< 

1 

0 

$* 


s i 

g 2 

% 


0 


\ s z 

®1 

S 2 


Y 

r- 2 

i ■ 

S 3 

CvJ 

tsD 

g l 

j 

1 

> 

Y 2 


10 0 


S> 1 +S 2 Y+g 3 Y 2 


S i +g 2 Y+S 3 Y 2 

.0 y " 1 0 1 


S 2 ^ ir ig 2Y 2 

=3 

„ ~-l ' 

« 2 Y + ®i + S 2 Y 

0 0 y “ 2 


e 3 + e 2 Y+s lY 2 

» ' 

^ 3 y"‘ 2 +S 2 +s iY 


.(4.67) 


Vq» ^1*^2 are com P lex -numbers. 

4.5 PERFORMANCE EVALUATION OF MISMATCHED RECEIVER . 

Proceeding as in the conventional receiver we can 


derive an expression for the performance index .and 


m 


IT\ 

V l K 2 Y 

*rh 

V i K 0 7 


.... (4.68) 


where ar0 NxM matrices whose elements are defined in 

equations ( 4.31 )» (4.34) and (4.35a). It can be verified that 
is real though V is complex. The performance of the 
mismatched receiver is numerically evaluated in chapter 5 using 
equation 4.68. 


4.6 ADAPTIVE MISMATCHED RECEIVER. 

The performance of the mismatched receiver is 
evidently non- inferior to that of the conventional receiver. 
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This is because even when the noise is predominantly white 
the mismatched receiver will adjust its coefficients such 
that it is nothing but a conventional receiver and hence it is 
optimum. As the noise becomes more and more coloured the ' 
mismatched filter coefficients change to combat the effect of 
noise (clutter) though not completly. For a conventional 
receiver the knowledge of the clutter characteristics is of no 
advantage. But a mismatched receiver filter coefficients can 
be adjusted with varying clutter characteristics to obtain 
bettor performance. Thus even in varying clutter environment 
the performance of mismatched receiver can be made superior 
compared to conventional receivor. 

The characterization of the clutter statistics by 

equation (3.6) is handy. The moment we estimate the parameters 

V ,Oa,V a the clutter characteristics are known. So the 

or’ r’ d’ od 

mismatched receiver coefficients can be computed. In a changing 
clutter environment it is enough to change the mismatched 
rocoiver coefficients only when the clutter characteristics 
change significantly. 0£in other words it is necessary to 
compute new set of filter coefficients only when the parameters 
change significantly. In this connection the parameters are 
quantized and the measured clutter characteristics is classified 
to that quantized parameter sot to which it is the nearest* 

For determining the quantization step for the parameters 
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V »V we specified the maximum degradation in 

or* od’ r* d ^ 

performance due to quantization as 0.1%. By trial arid 'error 
we arrived at the quantization step that will result in a 
performance degradation less than or equal to 0.1% .. The 
computer program output at some stage of the trial and error 
procedure is given in table 5.3 It is found that the parameters 
(V^d'Vod^ require (2, 8, 2, 8) quantization levels 
respectively. 


4.7 ESTIMATION OB THE CLUTTER PARAMETERS . 

The first stop in the adaptive receiver is the 
estimation of the parameters. The parameters a d' od arG 

to be estimated from the covariance f unction of noise, i.e. 
the covariance function of the received signal when the. target 
is absent. The covariance function of the noise 

v (t,u) * \ ( t,u)+No6(t-u) 
n q 

V (t-u) is given by (2.17) which wo reproduce here for 

i^c 

convenience 


K (t-u)= % K-(l+cosa'ce'' 1:rc ) e 3 c ‘ r ~ dT r(t) .... (4.71) 
n c. ** 


where a = 2V r x 2 - Tr 
b = 20^4^ 
c « 2itY d 

d= ic£x4--n x 


....(4.72) 


and R(t;) is given by 
equation (3*3) and (3*4) 
which is reproduced here 
for convenience. 
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R(%) 


; m-k 
■SL M 

k=o 


(1- 


6- I3& ) \ 

T > 


(4.73) 


where p(x k )= T-KT p if ..... (4.73 ) 

= 0 otherwise 

We observe that when x- IT p 


then R(x) = - 

2 

Then K n (NT p )= ^ K^( 2^^)(l+cosaxe" t)/r 


)e 


.... (4.74) 
2 

j CT-dT 


(4.75) 


where x=ITT p 


Prom the observed data the covariance function of the noise 
can be computed (W) and for delays x=ITT p ; N=C,1..M-1 
its value is given by eqn. 4.75* Using this K n (x) the 

C - 

parameters a,b,c,d are to be estimated. Since the unknown • 

parameters are exponentially related to (x) it is difficult 

c . 

to estimate tho parameters directly. 

Taking logarithm of (4.75) 

' bx^ 2 

log (NT p )= K 2 +log(M-N)+log(l+cosaxe~ )+jcT~dT ....(4.76) 

P 

For _ 

• *■■ 2 . .. 

1 °g el ( T p )=Z 2 +log( M-l )+log( 1+ c osaT p e P )+ 3 cT p -dT p ...(4.77) 

For W=2 


S3 


l°g ^ ( 2!Ep)=Z2+log(H~2)+log(l+cos2aTpe''^ T p)+3 cTp-dTp 

• * • •( 4.78) 

2 2 

Since cosaTp© P log(l+cosaT p e ^^P) can be expanded 

in power series. For simplicity we shall retain only 2 terms 
in tho expansion then (4.77) becomes 

• 2 2 - 5 

1 ° s E n c ^ T P ) =E 2 +l0 ^^ M “ 1 ) +x y Jr +3 cT p -dT p • * • • (4.79) 

* ■“ • *- 2 -■ cos2aT G “4bTp 

log K n (2T p )=K 2 +log(M-2)+cos2aT p e~ 41:)T p 

' + 2 j cTp-4dTp 

= K 2 +log(M-2)+(2x 2 -l)y 4 ~ +2 j cT p +4dT p 

. . . . .(4. 80 ) 

vm2 

where x= cosa5D p y=e P 

combining (4*79) and (4.80) we can get a multinomial equation 
in x,y eliminating the parameters c and d... However this 
procedure does not make use o:f the complete observed data and 
tho solution for tho multinomial equation is difficult to 
obtain. ■ 

However wo know from last section that the parameters 
^r'or'VJd) tak0 only ( 2 » 2 » 8 » 8 ) suanii 21 ** -levels. Hence 
(a f b,c,d) take only 256 combinations. Our job is to find 
which parameter set closely fits the observed data. From 
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(4*75) we observe that the parameter c can be evaluated, by 

ITT* Ya.+ ,o of " 

taking^ the real part of (t) by the imaginary part. Having 

' c 

found the value of c we are left with only 32 combinations. 

We have to decide to which parameter set of these 32 combinations 

'"iS 

the- measured data closest . For this we can adopt the 
following procedure : 

Let us find the covariance function of the received 
signal when the target is absent. Let the value of this for 
delays equal 0 , Tp... (M-l) Tp be arranged in a vector* 

Similarly using (4.75) for each possible parameter set covariance 
vector would be found* Then the problem becomes one of finding 
the nearest vector [12] to the observed vector from the 32 
vectors. By finding the distance wo can find the nearest vector* 

4.8 ADAPTIVE RADAR. RECEIVER SCHEME . 

We. shall summarize the above sections. The clutter 
parameters are estimated from the observed data as in the 
last section. The mismatched filter coefficients are procomputed 
and stored for each parameter sot. Hence once the parameters 
are knqwn the optimum mismatched Receiver filter coefficients 
are read. Since this schomo requires very little computational 
efforts it can adapt fast to thechanging clutter environment. 

4.9 DISCRETE RESOLUTION RECEIVER . 

The next receiver that is easy to implement that we 
consider is the discrete resolution receiver. The performance 
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vo have roplacod the dummy parameters m,k by m'k* in 
aquation (4.84) 

(4.84) can be rewritten as 
M-l 


f SZf(t-k»!D o -m*T p )f»(t~mT ‘ dt 

J TQrrO ^ 


nt=o 
M-l N-l 


+ J 2T SZl i, (t-M o -“ I p) f ( 1; - k 'V m,T p)6L«=0 .... (4.85) 

•m=o k=o 

For a rectangular pulse defined by (3.1) 
f ( t-kT 0 -mT p )f* ( t-k’ T 0 -m » T p )= .... (4.86) 

4. mT 

J f r (t-ma^)f(t-k"T o -m*ai p )e“ 3 Vdt= £ 2 -^-- ;P 6(k* )6(m-m* ) 

..,.(4.87) 

and J" f ( t~kT 0 -mTp )f ( t-k * $ 0 -m ’ Tp )-6 ( k~k 1 ^(m-m* )....(4.88) 

substituting (4.86),(4.87)and (4.88)in (4.85) 

M-l u M-l 3ST-X 

^^jvTpm 6(k , ) + jg~ Tl^Oc-k 1 )6(m-m« ) 

n»o m — 0 k=o 

» J^(t)f(t-k , T 0 -m‘I p )dt .... (4.89) 

♦> 

For various values of k* »m’ (4.89) becomes 

g* 0 + 1 * jV r (*)f( t ) dt 

:... (4.90) 

s 20 = 


g* ( t )f ( t-2T p )dt 

(t )f ( t-IPlTp )dt 

^01 43 s rS* (t )f(t-! 0 )dt 

# - f|(f )f(t-2T 0 )dt and so on 


%-10 
«02 


Since g(t)= J<Hf )e*“ 32wft dt ....(4.91) 

Substituting G(f) from equation (3.58) we get 
- \ 

Jg» ( t )f ( t-kl 0 -mlp )dt= JJ e “3 2%ft df f ^ ) dt 

.... (4*92) 


= iV) 6 ' KfcVmTp)wp(f) 


93) 


-g p | ^ ^ ( f )e j *’ T p( W+V D ) exp(-3(kT 0 +mT p )w) 

• • • ♦ ( 4 • 94 ) 


= ^ X ei<!VD p f {f+f D )F(f)9 

ifc 0 i? • — 


j(m-£ T p +kT Q )w 


•< “n^ 

*0 


where F( f )= J"f( t)e’"^ 27lft dt 


df ....(4.95) 


(4.96) 


— 


So the main task in evaluating the values of is to find 
the integrals of the form 
- r jwT n P 

df *...(4.97) 


j /(f+f d )P(f)e‘ 




«° s n (f) 


■ r* t 

1*, s» o W' 


where P=0,+1 » ♦ ...MN~1 
2H-MT-1 


df 


.... ( 4 . 98 ) 


where S (f) is defined in equation.. 
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Closed, form solutions for the integrals given by 
(4*98) donot exist* However we can evaluate these integrals 
numerically. We shall consider the numerical procedure in the 
next chapter. 


4.10 PERFORMANCE EVALUATION OF DISCRETE RESOLUTION RECEIVER * 

The receiver correlation waveform corresponding to 
discrete resolution receiver is given by 


M-l 


M-l M-l 


s(t)« JjeZ f ( t-mTp )e~^ w D^’+ 51 5” f(t-kT -mT Jg^ •♦..(4.99) 
r m=o k=o r ELK ’ 


m=eo 


let Tp=NT o * The expression for the performance index 
^ can bo derived as follows: 

O' 


Now L 


J r(t)g^(t) dt 

T T+T (M-l)T p +(N-l)T 0 +T. 

r(t)^(t)dt+ J~ r(t)g*(t)dt +*...' J' 


0 


“ T. 


(M-l )Tp+( N-l )T ( 


r(t)^(t) dt 


T 


(M-l)T +(N-1)T +T * 

1 . 0 . . 6 M-lN-3 

.fT 


r, [e^D^+gnn] f . v . e M-lN-l 

S r(t) p— J dt+... J - *<*> dt 


(M-l)T p +(N-l)T 0 

(4.99a) 


» • « # 



V 
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*fc 

Again e D‘ can be assumed to be constant over a period of T 
Then 


^ C 0 °1* •• #C MN-l] 


U 


0.. 


u. 


u. 


MN-1 


• « • >( 4< 100 ) 


raT^+nT +T 
o 


where C. 


■ ■“-p 

i=/ 


r(t )dt 


mT D +nT. 
P o 


to = imod M 1= 3 mod M 


i~K . .1«1 

n = -g— n= -jf— 


and TJ = 


•‘Tw 


mr» 


YT 


. ’.J^ 2 - n * T c^ t cu 

% l» ' to. 4 

2 /H 1 3=U ,(T E[C <t!P C3 U 4 u* T VU where Y=E[C* T C3 
Now V(i,j) = E[Cj* C.j3 

mTp+kT Q +T nTp+i-T Q +T 

« J J B[r*(t)r(-u)]dtdu *■ . 

mT p +kT 0 nTp+ T Q s 

= (T [K (t,u)+No6(t-u)4Blfje B ('fc)f e (u)e"3 w D(' t - u )3 ( itdu 
mTp+kT Q +T »Tp+ T Q +T 

Now J' j' No6(t-u)dtdu = NoT6(m-n)6(k~£) 

m!Ep+kT j0 nTp+ T Q 


sec. 
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mTp+kTo+T aT p +|T 0 +T 

i J f (t)f s {u)e^ w D (t ~^/fB^7dtdU 

mT P +kT o nT p + T o 

= TeL|B( 36(^.)6(k)e'“ ;3w D( n “ n ) T p 
J"f K n c ^ t,u ^ tdu ^ dT p -nT p ), JJ R( t -u )d tdu 

• -^(5^ T p ) f T6(k~£ )( ^p3L) 

7(1, j)= | aCk-OK^CSPS!^) [ Ji “' [ ni “ n) l+IToT6( m -n)6(k-.i!.) 

+ T£lBl3fe(£)6(k)e^ w D (m " n)T P- 
••■ A- jQjUL ■ .....(4.101) 

Y tJ 

•where X,Y are MxM matrices whose elements are defined by 
X (i r j) = !I?6( J g)(k)e'* ;3 V m ~ n)T P 

Y (1,3) » 7(i,j) - X(i,j) 

UBing' equation (4.101} the performance can he 
evaluated. We shall consider this in the next chapter. 



In this chapter the results of the numerical 
evaluation of the performance of the receivers considered in 
the previous chapters are presented. The listing of the 
programmes executed on DEC— 10 computer are included at the end 
of the thesis. Along with the presentation of relevant 
programme outputs .assumptions made in the performance evaluation 
are briefly discussed. 

5.1 NUMERICAL EVALUATION of THE PERFORMANCE 0F r _THE. OPTIMUM 
RECEIVER . 

The expression for the performance index for the 
optimum receiver is given in (3.62) which is reproduced here 
for convenience 


06 p 

A w- f 

o = ^ J No4-S (f 5 f 


.... (5.1) 


where F (f-f D )= ${72 t-.Z'’’ T~ ' ••••(5.2) 

B jj • 1 'l * 7i * 1=0 


We evaluated (5.1) using numerical integration. The step size 
had been chosen so as to sample all the frequency components 
adequately. Since, about 90* ofthe area of the Sine 2 curve 
is contained between - ^/the limits of the integration in (5.1) 
|re kept at - -f and f and the integer H is chosen by the 
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following consideration. We know that the performance index 
corresponding to the case when clutter is absent is given by 

sl * ¥e com P utod ^ 0 by numerical integration of (5.1) with 
i value of N selected so as to make the difference between A 
md less than 2.5/^. 

For the computations we assumed the pulsewidth, 
Lnterpulse period to bo 1 pSec. , seconds respectively. 

?or 97% accuracy in the computation I was found to be two for 
* kits/scan and the number of numerical integration intervals, 
*idth\of the interval were 8x10^, 4 respectively. 

>•2 RESULTS A HD OBSERVATION . 

The performance indices for the conventional, 
aismatched, white-noise optimum, optimum receivers are computed 
ising equations (4.35), (4.58) and (5.1) respectively. 

The number of pulses in the transmitted waveform 
5 hat is reflected by the target over one scan period 
[no of hits/scan) is taken to be four for comparing the 
performance of the conventional, mismatched receiver with the 
optimum receiver. The results of the computation are 
tabulated in Table 5.1. 
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In tablo 5.1 tho notations used have the following meaning. 


SYMBOL 

SIGNAL 

WHITE 

CLUTTER 

TDOPPLER 

CDQPPLER 

IDEAL 


MEANING 

TARGET SIGNAL power/ Sample 

Spectral density of white noise 

Cluttor powor/sample 

Target doppler frequency 

Doppler freq. due to clutter 

Aj-THE PERFORMANCE INDEX FOR 
OPTIMUM RECEIVER IN THE Absence of 
clutter. 


MATCHED for conventional receiver 

MISMATCHED for mismatched receiver. 


From tablo 5*1 wo observe tho following : 

1) ^ or tho white-noise optimum (denoted as ideal) receiver 
should be 40 for the given signal and noise parameters. 

But A ' , the value of ^ computed using numerical procedure is 
38.97. Hence the numerical procedure has about 2.53£ 


inaccuracy. 


2) For tho same signal and white noise powers. She performance 
of optimum receiver is superior to be other two receivers for 
all clutter powers and for all target doppler shifts, as it 
should be . 


5) The superiority of the optimum Rx increases with increasing 
flutter power. This may be verified by computing Q > 

rphe result of such an analysis is given in Table 5.2 


or some typical values from table-5.1* 
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TABLE 5.2 


COMPARISON OP SUB OPTIMUM RECEIVERS PEREQRMANCE INDICES 
WITH THE OPTIMUM RX. 


TARGET DOPPLER = 100 Hz 
* 

Target Signal Power= 10 

White noise No = 1 
Spectral density 

B = 6261. 8793 

B+D= 9044. 9368 

A =111. 9096 





mismatched? % improvement in in 
xoptimum Rx com pared tc 
xConvent- /^ematcnea - 

__.Jiio.nal Rac_ t 


1.8 


1.4 


CLUTTER 

POWER 


0.1 

36.93 

1.0 

26.17 

10,0 

8.95 

100.0 

2. 06 

1000.0 

0.34 

10000.0 

0.04 


36.26 

36.39 

19.69 

21.74 

3.53 

5.17 

0.38 

0.62 

0.03 

0.06 

0.003 

0.006 


32.9 

20.3 

153.5 

73.1 

442.1 

232.2 

1033.0 

466.0 

1233.3 

566.0 


4) At large dopplers the difference in performance of the 
above receivers is small. 
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We obtained the receiver operating characteristics 
viz the plot of probability of detection. Vs' prpbabilit^r.. ."p 

of False alarm and are given in Fig. 5. 1-5 . 6 for some 

representative values. 

The computational burden involved in evaluating the 
performance of the optimum receiver for the case of 16 hits 
per scan is excessive. Therefore, for the case of 16 hits per 
scan. Wo limit our comparison of performance of the above 
structured receivers’ with that of the receiver that "would be 
optimum in the absence of clutter. The results of the 
computa,tion are tabulated in table 5.3. From the table 5.3 
we observe the following : 

1 ) When the relative doppler between the clutter and the 
target increases the performance improves. 

2 ) We observe that increasing target power improves the 
performance . 

3) When the noise is predominantly white both the 
receivers have the same performance. Otherwise mismatched 
RX is superior to the conventional R X. 

4) As'the clutter power increases the superiority of the 
mismatched Ex over the conventional Rx improves* 

ROC’s for some typical values of A are shown in 


Fig. 5.7-5.10 
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We have not attempted to compute the performance 
of the discrete resolution receiver due to the relatively 
large computational time requirements. 
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CHAFPi8-i . 

CONCLUSION AND SUGGESTIONS FOR FURTHER WORK 
6,1 CONCLUSION. 

In this thesis we studied some receivers for a 
pulsed radar to detect a Swerling 1 target in the presence of 

rh«c!utt«v t-viociel . . 

a random dipole modelled clutter^incorporated xs representative 
of chaff and vegetation clutter characterized "by a spectrum 
with a secular component arising from the rotationof the 
scatterere. For the pulsewidths and pulse repetition frequencies 
normally used in a pulsed’ radar , the aircraft can he modelled 
as a slowly fluctuating point target and we assumed a 
Swerling 1 target. 

For the above target detection problem we attempted 
to solve the associated integral .equation. Although we succeeded 
in finding a solution involving complex error functions, we 
found it too unwieldy to be of any practical interest, be 
then derived an approximate solution under the realistic 
assumption t ha t the clutter correlation function can bo treated 
as a constant over intervals of one pulsewidth of the trans- 
mitted pulse train. However this solution only gives us the 
Fourier transform of the time signal whose cross correlation 
with the received signal the optimum receiver would evaluate 
before thresholding and making a decision about the presence 
or absence of the target. As an analytic expression for the 
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inverse Fourier transform which would have in an explicit form, 
the dependance of the required time signal on the clutter 
parameters, escaped our efforts we did not consider any further 
the implementation aspects of the optimum receiver. However, 
to allow for a comparative evaluation of the performance of 
other receivers with that of the optimum receiver, we obtained 
an expression for the optimum receiver performance. 

In view of the above observations about the optimum 
receiver we next considered three relatively easily implementabl 
structured- optimum receivers. They are the conventional 
receiver, the mismatched receiver and the discrete resolution 
receiver. An adaptive version of the mismatched reccivor was. 
proposed for changing cluttor environment. In this connection 
the parameters ( a ^* G T *^ 0T > v oa ) wore quantized and for onsuring 
the degradation in performance duo to quantization to be below 
0.1 X it was found that (8, 2,. 2, 8) levels were sufficient. 

We concluded that the parameters can:., be estimated by computing 
the clutter spectrum and then using tho nearest neighbourhood 
method. , 

The performances of the mismatched and conventional 
receivers were evaluated for various combinations of signal , 
white noise cluttor power and for different clutter and 
target dopplers for 16 hits/scan. Their performances wer© 



compared with the white-noise optimum receiver in the 
absence of clutter. The mismatched receiver performs^ better 
than the conventional receiver at large clutter power. Mie'n 
the noise is predominantly white their performances are the 
same. Their performances were .compared with the optimum 
receiver in the presence of clutter at 4 hits/scan. Their 
performance is * satisfactory’ at large target dopplers. The 

t 

* cii&cv'ate YQSoIiaWovs 

performance of the receiver could not bo evaluated 

because of excessive computational requirements. 

6.2 SU&GKBSTION . FOR FURTHER WORE. 

As an extension of the work done in this thesis 
a comparative study of the receivers proposed in this thesis 
and that considered by Haykins et al [ *9 3 may be undertaken.. 
Haykins et al [ fS 3 uses the random dipole clutter model for 
designing an MTI Radar receiver with constant False Alarm rate. 
If the nth order FIR filter used in this receiver is designed 
to maximize the probability of detection performance then such 
a comparison is feasible. 

Further, if the clutter parameters are quantized then 
the solution to the integral equation associated with tho 
optimum receiver may be numerically evaluated to store the 
corresponding optimum filter coefficients for each possible 
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combination, of the quantized parameter values. Then in a 

changing 'clutter environment as soon as the clutter parameters 
are estimated the optimum filter coefficients may he used to 
process the received signal. 

In stead of using the random dipole model for clutter 
directly as we did in this thesis and as suggested above one 
may fit an Nth order ARMA model to the clutter. Hie value of 
N for representing the clutter covariance within reasonable 
accuracy may be found to implore the implementability and 
adaptability of the associated whitening filter. 
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* u * ( :> • 1 4*54O*.0*O.O33)**2 
rc*2«O*3«l4*Q*t*540 f 0 
B^ns^.O*? 3. 14*540. 0*0. 044)**? 

ACCEPT */;n^2,N3 J rN4iN5,Nb,N7,N8,N9 <- 
M 1 • N ? GIVES THE RANGE OF VAL ,: t.S n v?h 

Wllf CM the PARAMETER a will bf: vapifj 
— ii 3 . w 4 ■■•Range for d 
•’S,N6*RAN(JE FOR A 
M7 JhR-RANGE FUR C 
^ m || - 
nu 10 I.*N1 ,N 2 

4L =L-l 
RaRB*AL/8.0 
no 10 iTaM 3 , N 4 

od = CRBotaJ' 1/4 .0 )+ti 
no 10 KsN5 
AK»K-1 f . 

no 10 n*n?*n8 

6«JSN-1 

i'=AA*.'.f f /R.O 
v- =1 ,y/54o.o0 

ti i Pi=tj*TP*TP 
r. T oisfin*Tp*TP 

Ai dba*TP 

opps-r+Tp 

M = M+ 1 

nu 10 T*1 #NuP 

*0*( »S1)T*B1P1) 

5!ii;»l4jitKKfiaSi))*co S ux«»wi 

rSSSBfl'sftUfiiS^mofe of mis roKWiEc 


f u=g-o. 

ssttiO. 

n r=i,i6 . 

wyvWfc i-JMt - .f«i «#• 


mm 



... "V:\ vjp ? C I iu}) ) B 1G~P C I .# J ) 

TF CS* AM, , GT . p C I r j ) ) SM ALl = p (. T . J ) 

(SufiTTSME 

nrd- ,A!1=0, 044*540. U*AJJ<B#0 
M'' : AL=0 .033*540. 0*AL/8.G 
•i ■■'. r v =o. 033*540. 0*AK/8.0 
■ ^0.28 
v,!si!T J 
VI -<> 

Y2 = 1 k 

r, r.. T * f *l,X2 

P,t T.*T » , L , i 2 

nuT.p’ 1 o 1 , l i' ( 1 »•■-)» V = 1 > 8 ) 
(.‘":i T I , !jr 


M)3 
1 03 u 


fln ! : r \ \ I (\ I . I 0 ? . 1 *'4 f SO ) 
r.- r’^TiiT 1 (>?(> .oTG^aO , finFREO » BPFPrlO , 

,?x! '! Z S IGH A n = ' VF b 1 2 V 1 t-Rl PEC = ’' ?Fb . 2 . 1 X / ' OPF 

ouTiVt io30*STGMAR,PHFRF.O#l>rFFE9 
-1 P ?SlMRR= , ,F?“a^-OW*EB. s lF5.?.7X.'PPFR E 4=>,f5. 


OPFRE0=' r F5.2) 


2 ) 


r.(, Tn 7 

PHl.iT luS,fHGi4AD,i>IGMAR,DPFREU 
FljWHM'ClX, ' VARYING DHFREg.HFANf* , 

2X. 'SIGMADs' ,F5.2,2X, 'STGMAP::' ,F5,2,IX, 

GO Tn 7 

OHTHT ,STGMAD, SIGMAR,DRFREO 
FURtlM'f IX, ' VARYING DUP.FREU. MEAN', 2X, 
*SIGMAD- f F5 .2 » 2X» SIGMA K®' ,F5.2,2X r 


'DPFREOs' /E5.2) 
# ORFRE0=s',F5.2) 


PHT.iT 1U6 

r<.Ri.AT(5X, 
f’AEl. GRAPH 
STOP 
FWO 


'NORMALIZED delay 


2GX, 'CORRELATION') 


FlPbT DIMENSION OF Y-UIMFNSinN OF X 
SECOND DTU. OF ¥=HO.OF GRAPHS PEP PAGE 
TIi THE INPUT FILE NAMED FOPLPT.DAT THE 
lsf LIME SHOULD HAVE THE Min, VALUE OF 
X, MAX. VALUE. 2nd LINE SHOULD HAVE YMin, 

THF A i’aST BUT ONE LINE SHOULD HAVE THE 
TITLE OF THE GRAPH IF ANY. THE LAST 
SHOULD HAVE TITLE FOR XAXIS.1Q BLANKS, 

TITLE for YAXIS 

DIMENSION XCH)»Y(U,4),AUS8). bl 51 ) 

■npFN ( UMIT»6 » DEV ICE® ' DSK ' ,FILE= P AMES .-jAT ) 

* pFn (tJNIT=22,DE VICE® ' DSK *, FILE® 'FuRlPX.DAT' 
‘ )AC147) ,A(148) 

)A(145),A(146) 


( Y ( I , J ) , J=l,4) 


f? i£ A \j C 2 *?> f 
r 1&L(22,». 
no 44 T=1 ,11 
READ ( 22 ,*) X ( I ) 
continue 
READ(2?,10)(A(I),I»1,144) 

FORMAT C 72 Al ) 

A { 1 4R )®0, 0 

CALL USPLX ( X , Y ,11, 4, 1,11, A , I PAG 4 , LF-O 

ffWTINUE 

CLOSE (UNI T»b , DEVICE® ' DSK ' .FILE® 'PANES . DAT ' 1 
CALL UERTST 

Call usmnmx 

STOP 

END 



um 


, ££ R£r c 


